import numpy as np


# 路面磨耗指数（PWI）
# PWI - 100 * a0*WR**a1
# WR = 100 * (MPDc - min(MPDl, MPDr))/MPDc

# WR 路面磨耗率（%）；
# a0 模型参数，采用 1.696；
# a1 模型参数，采用 0.785；
# MPD 路面构造深度（mm）；
# MPDc 路面构造深度基准值，采用无磨损的车道中线路面构造深度（mm）；
# MPDl 左轮迹带的路面构造深度（mm）；
# MPDr 右轮迹带的路面构造深度（mm）。


def WR():
    MPDc = np.random.randint(15, 18)
    MPDl = np.random.randint(13, 15)
    MPDr = np.random.randint(13, 15)
    return 100 * (MPDc - min(MPDl, MPDr)) / MPDc


def PWI():
    a0 = 1.696
    a1 = 0.785
    wr = WR()
    return 100 - a0 * wr ** a1


if __name__ == '__main__':
    pwi = PWI()
    print(pwi)
